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Abstract 



Several cosmologically distant astrophysical sources may produce 
high-energy cosmic neutrinos (E > 10 6 GeV) of all flavors above the 
atmospheric neutrino background. We study the effects of vacuum neu- 
trino mixing in three flavor framework on this cosmic neutrino flux. 
We also consider the effects of possible mixing between the three active 
neutrinos and the (fourth) sterile neutrino with or without Big-Bang 
nucleosynthesis constraints and estimate the resulting final high-energy 
cosmic neutrino flux ratios on earth compatible with currently existing 
different neutrino oscillation hints in a model independent way. Further, 
we discuss the case where the intrinsic cosmic neutrino flux does not 
have the standard ratio. 

PACS number(s): 14.60.Pq, 98.70.S, 29.40.K, 98.54.Cm 
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I. INTRODUCTION 



High-energy neutrino (E > 10 6 GeV) astroparticle physics is now a rapidly developing 
field impelled by the need for improved flux estimates as well as a good understanding 
of detector capabilities for all neutrino flavors, particularly in light of recently growing 
experimental support for flavor oscillations 

Currently envisaged astrophysical sources of high-energy cosmic neutrinos include, for 
instance, Active Galactic Nuclei (AGN) and Gamma Ray Burst fireballs (GRB) 0. Produc- 
tion of high-energy cosmic neutrinos other than the AGNs and GRBs may also be possible 

!• 

High-energy neutrino production in cosmologically distant astrophysical systems such as 
AGNs and GRBs follow mainly from the production and decay of relevant unstable hadrons. 
These unstable hadrons may be produced mainly when the accelerated protons in these 
environments interact with the ambient photon field (jry) and/or protons (pp) present there. 
The electron and muon neutrinos (and corresponding anti neutrinos) are mainly produced in 
the decay chain of charged pions whereas the tau neutrinos (and anti neutrinos) are mainly 
produced in the decay chain of charmed mesons in the same collisions at a suppressed level 

I- 

Previously the effects of vacuum neutrino mixing on the intrinsic ratios of high-energy 
cosmic neutrinos, in the context of three flavors, are briefly considered in ||. Here we update 
the description and consider three and four neutrino schemes with the most up-to-date 
constraints from the terrestrial, solar and atmospheric neutrino experiments. We also discuss 
the case where the intrinsic high-energy cosmic neutrino flux has nonstandard ratio which 
may not be obtained in charged pion decays. 

The present study is particularly useful as the new underice/ water Cerenkov light neu- 
trino telescopes namely AMANDA and Baikal (also the ANTARES and NESTOR) as well 
as the large (horizontal) shower array detectors will have energy, angle and possibly fla- 
vor resolutions for high-energy neutrinos originating at cosmological distances J7|. Several 
discussions are now available pointing towards the possibility of flavor identification for 
high-energy cosmic neutrinos ||. 

The plan of this paper is as follows. In section II, we discuss the effects of vacuum 
neutrino mixing in three as well as four flavor schemes and numerically estimate the final 
(downward going) ratios of flux for high-energy cosmic neutrinos using the ranges of neutrino 
mixing parameters implied by recent searches for neutrino oscillations. In section III, we 
discuss the consequences of a hypothesis in which the intrinsic cosmic neutrino flux does not 
have the standard ratio [F°(u e ): F Q {v^)\ F°(v T ) = 1: 2: 0]. In section IV, we summarize 
our results. 
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II. EFFECTS OF NEUTRINO MIXING ON HIGH-ENERGY COSMIC 

NEUTRINO FLUX 

In pj and in pp collisions, typically one obtains the following ratio of intrinsic high-energy 
cosmic neutrinos flux: F (u e ): F (u/j): F {y T ) = 1: 2: < 1CT 5 . For simplicity, here we take 
these ratios as F°(v e ): F°{u^)\ F°(u T ) = 1: 2: 0. We also discuss the effects of vacuum 
neutrino mixing by varying first two of these ratios from their above quoted values as this 
might be the case under some circumstances ||. We count both neutrinos and anti neutrinos 
in the symbol for neutrinos. 

We consider an order of magnitude energy interval essentially around 10 6 GeV since the 
currently available models for high-energy cosmic neutrinos give neutrino flux above the 
relevant atmospheric neutrino background specifically from AGN cores within this energy 
interval. Also in this energy range the flavor identification for high-energy cosmic neutrinos 
may be conceivable in new km 2 surface area neutrino telescopes ||. We have checked that 
for high-energy neutrinos originating from cosmologically distant sources, the change in the 
flavor composition of the high-energy cosmic neutrinos due to vacuum mixing is essentially 
energy independent for the entire energy range relevant for observations as the energy effects 
are averaged out in the relevant neutrino flavor oscillation probability expressions. 

It has been pointed out that there are no matter effects on vacuum neutrino oscillations 
for relevant mass squared difference values [(9(10~ 10 ) < Am 2 /eV 2 < 0(1)] for high-energy 
cosmic neutrinos scattering over the matter particles in the vicinity of sources [Hj. Never- 
theless, for high-energy cosmic neutrinos scattering over the relic neutrinos, in halos or other 
wise, during their propagation, there may be relatively small (at most, of the order of few 
percent) matter enhanced flavor oscillation effects under the assumption of rather strong 
CP asymmetric neutrino background |Trj . However, given the current status of high-energy 



cosmic neutrino flux measurements, vacuum flavor oscillations remain the dominant effect 
and therefore from now on we consider here the effects of vacuum flavor oscillations only. 



For some other possible mixing effects, see fl2fl . The typical distance to these astrophysical 
sources is taken as L ~ 10 2 Mpc (where 1 pc ~ 3 ■ 10 18 cm). It is relevant here to mention 
that our following analysis is not necessarily restricted to high-energy cosmic neutrinos orig- 
inating only from AGNs or GRBs as the above mentioned ratios of intrinsic neutrino flux 
can in principle be conceivable from some other cosmologically distant astrophysical sources 
as well 0. For simplicity, we assume the absence of relatively dense objects between the 
cosmologically distant high-energy neutrino sources and the prospective neutrino telescopes 
so as not to change significantly the vacuum oscillation pattern. 
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A. Three neutrino scheme 



l(T 10 eV 2 and Am 2 arger 
jAm^l < |Am§ 2 | < [Amf-J where Am^ 



It has been known in the two flavor framework that the solar and the atmospheric 
neutrino problems are accounted for by neutrino oscillations with Am 2 mallcr ~ lCT 5 eV 2 or 

10 _2 ' 5 eV 2 , respectively. Without loss of generality we assume that 
mf — m 2 (i ^ j; i,j = 1, 2, 3). If both the solar 
and the atmospheric neutrino problems are to be solved by energy dependent solutions, we 
have to have Am^ ~ Am|, and Am| 2 ~ Am^ m , i.e., we have mass hierarchy in this case. 
Here, the subscripts and atm stand for the mixing angles for the solar and atmospheric 
neutrino oscillations, respectively. 

In the three flavor framework, the flavor and mass eigenstates are related by the MNS 



matrix U 13 
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In the present convention for the mass pattern, the disappearance probability for the CHOOZ 



experiment [14] is given by 



l-4|f/ e3 | 2 (l 



(2) 



and combining the negative result by the CHOOZ experiment [13] with the atmospheric 
neutrino data, we have (for quantitative discussions, see WU): 



u e ,v < i. 



(3) 



As was mentioned in the Introduction, the matter effect is irrelevant in our discussions, 
so let us consider the vacuum oscillation probability in the case where the oscillation length 
is very small as compared to the distance between the source and the detector [hereafter 
referred to as far distance (approximation)]. In vacuum, it is well known that the flavor 
oscillation probability is given by 



P{y a -> vp\ L) = S, 



In the limit L — ► oo, we have 



■iAE jk L 



(4) 



P(u a -> Up] L = OO) = S af3 - U aj U t3jU ak Up k = \ U aj\ 2 \Upj\ 



(5) 



where we have averaged over rapid oscillations. Thus, we can represent the oscillation 
probability as a symmetric matrix P and P can be written as a product of a matrix A: 
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Now, the cosmic neutrino flux in the far distance can be expressed as a product of P and 
the intrinsic flux F°(u a )(a = e, /x, r): 



/ F(u e ) N 



F{u T ) 



( po( v \ \ 



P 



J 



F°(u„ 



1 F°(u P ) ^ 



AA 



T 



J 



F>») 
F°(u T ) 



(8) 



Throughout this section we assume the standard ratio of the intrinsic cosmic neutrino flux, 
i.e., F°{u e ): F°{u T ) = 1: 2: 0, so that, we get 
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where we have used the unitarity condition, i.e., Y^j \U a j\ 2 — 1- When \U e3 \ 2 <C 1 and 
\Ufjsl ~ \U t3 \, we have ||£7 W -| 2 — |^rj| 2 | "C l(j = 1,2,3), so the second term in Eq. @ is 
small. Hence with the constraints of the solar and atmospheric neutrino and the reactor 
data, we obtain 
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where we have used the unitarity condition again. Therefore, we conclude that the ratio of 
the cosmic high-energy neutrino fluxes in the far distance is 1: 1: 1, irrespective of which 
solar solution is chosen. This is because the MNS matrix elements satisfy \U e3 \ 2 <C 1, 
{U^l — \U T s\ and F (u e ): F°{i/^): F°{y T ) = 1: 2: 0. This is the feature that is expected to 
be observed at the new km 2 surface area neutrino telescopes in the case of the three neutrino 
oscillation scheme. 

Using the allowed region for the MNS matrix elements given in WE |, we have evaluated the 
final ratio of the cosmic neutrino flux numerically. To plot the ratio of the three neutrino 
flavors, we introduce the triangle representation, which is introduced by Fogli, Lisi, and 
Scioscia [16] in a different context. In Fig. 1, a unit regular triangle is drawn and the 
position of the point gives the ratio of the final high-energy neutrino flux with F a = F(u a ). 
One reason that we adopt this representation is because currently we do not know the 
precise total cosmic neutrino flux, while we may observe the ratio of flux of different cosmic 
neutrino flavors experimentally to a certain extent. Using this triangle representation, the 
allowed region in the three flavor framework with the constraints from the terrestrial, solar 
and atmospheric neutrino data, is given in Figs. 2a and 2b. The allowed region is a small 
area around the midpoint F(u e ) = F{v^) = F(u T ) = 1/3 and the small deviation from the 
midpoint indicates the smallness of \U e3 \ and ||i7 w -| 2 — |t^ T j| 2 |- 



B. Four neutrino scheme 

Let us now consider the case with three active and one sterile neutrino. Schemes with 
sterile neutrinos have been studied by a number of authors [I7|-^]. Here, we are interested 



arc 



in the four neutrino scheme in which the solar, atmospheric and LSND experiment [2T 
all explained by neutrino oscillations. By generalizing the discussion on Big-Bang Nucle- 
osynthesis (BBN) from the two neutrino scheme |22| to four neutrino case, it has been shown 
|19| , |20[| that the neutrino mixing angles are strongly constrained not only by the reactor data 
|T8| but also by BBN if one demands that the number N v of effective neutrinos be less than 
four. In this case, the 4x4 MNS matrix splits approximately into two 2x2 block diagonal 
matrices. 



On the other hand, some authors have given conservative estimate for N v |23j and if their 
estimate is correct then we no longer have strong constraints on the neutrino mixing angles. 
In this section, we discuss the final ratio of the cosmic neutrino flux components with and 
without BBN constraints separately. Throughout this section we assume the following mass 
pattern 

ml < ml <C ml < m 2 , (11) 

with corresponding Am^, Am 2 3 and Am 2 2 stand for the mass squared differences suggested 
by the solar, atmospheric neutrino deficits and the LSND experiment, respectively. 
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1. Four neutrino scheme with BBN constraints 



Here, we adopt the notation in |L9| for the 4x4 MNS matrix 
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where Cy = cos Sy = sin^- and 



(12) 



(13) 



Rjk(O) = exp (iT jk 9) , 
is a 4 x 4 orthogonal matrix with 

(Tjk)i m = i (SjiSkm — 5jmb~ke) , (14) 
and 2A 3 = diag(l, —1, 0, 0), 2y / 3A 8 = diag(l, 1, —2, 0), 2y / 6Ai 5 = diag(l, 1, 1, —3) are diago- 



nal elements of the SU(4) generators. From the constraint of the reactor data of Bugey [24 
which strongly constrain the disappearance probability P(u e —>■ u e ) for the entire region of 
mass squared difference implied by the data of LSND |2l| and E776 f2"5H , we have Hl8|-p0 



U e3 \ 2 +\U e4 \ 2 < 1 



(15) 



On the other hand, from the BBN constraint that sterile neutrino be not in thermal equi- 
librium, we get 



\U s3 \ 2 + \U si \ 2 < 1. 



(16) 



In this case, the oscillation probabilities in the far distance are given essentially by the 
following two neutrino flavor formulae: 
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where, we have substituted 6*13 = 9 14 = #23 = #34 = in Eq. (|T|) and #12, #24 correspond 
to the mixing angles of the solar and atmospheric neutrino oscillations, respectively. It has 
been known |26j that only the Small Mixing Angle (SMA) MSW solution is allowed in the 
v e *-* v s solar oscillation scheme, so in the present case, we take \9\2\ <C 1. Thus, we have 
the following ratio of the final cosmic high-energy neutrino flux: 
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2 - sin 2 29 atm 
sin2 2 #atm 

V ° J 

where, we have taken # 2 4 = $atm which satisfies [27] 

0.88 < sin 2 2^ atm < 1. 

In this case, therefore, we again have F{u e ): Fiy^: F(u T ) 

We have also numerically obtained the allowed region by letting 9i 2 , #13, #14, #23, #34 and 
#24 to vary within the constraints obtained from the reactor, solar and atmospheric neutrino 
data. In the four neutrino scheme, the total flux of active neutrino is in general not the 
same as that at the production point, and in principle we have to use a unit tetrahedron 
to express the ratio of the four neutrino flux. In practice, however, we may not observe the 
cosmic sterile neutrino events nor do we know the precise total flux of the cosmic high-energy 
neutrinos, so it is useful to normalize the flux of each active neutrino by the total flux of 
active neutrinos: 



(19) 
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(20) 



After redefining the flux, we can plot the ratio of each active neutrino with the same triangle 
graph as in the three neutrino case, and the result is depicted in Fig. 3. The allowed region 
is again small and the reason that the region lies horizontally is because of possible deviation 
of 9 2 4 = #atm from [ see Ec i- AH)]- 
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2. Four neutrino scheme without BBN constraints 



In this subsection, we discuss what happens to the ratio of the final cosmic high-energy 
neutrino flux if we lift BBN constraints. Without BBN constraints, the only conditions 
we have to take into account are the solar and atmospheric neutrino deficit data and the 
results of experiment of LSND (the appearance experiment of — > u e ) pi |, Bugey (the 
disappearance experiment of v e — > i> e ) and CDHSW (the disappearance experiment of 
Vn) P5|. For the range of the Am 2 suggested by the LSND data, which is given 



by 0.2< AmL SND /eV 2 ^ 2 when combined with the data of Bugey [f2|] an d E776 the 
constraint by the Bugey data is very stringent and 



U e3 \ 2 + \U e A 2 < 1Q- 2 



(21) 



has to be satisfied |T8|- pO|1 . Therefore, for simplicity, we take U e3 = t/ e4 = 0, in the following 
discussion. 

The analysis of the solar neutrino data in the four neutrino scheme with ansatz U e3 = 
U e 4 = has been done recently by Giunti, Gonzalez- Garcia and Pena-Garay |[29|| . They 
have shown that the scheme is reduced to the two neutrino framework in which only one 
free parameter c s = \U s i\ 2 + \U S 2\ 2 appears in the analysis^ Their conclusion is that the 
SMA (MSW) solution exists for the entire region of < c s < 1, while the Large Mixing 
Angle (LMA) and Vacuum Oscillation (VO) solutions survive only for < c s < 0.2 and 
< c s < 0.4, respectively. 

The analysis of the atmospheric neutrino data in the four neutrino framework has been 
done by one of the authors JJ(J more recently again with ansatz U e3 = U e A = 0. The 
conclusion of |3D| is that the region — n/2 < 6*34 < it/ 2 and < #23 ~ tt/6 as well as 
< 5\ < 7r is consistent with the Superkamiokande atmospheric neutrino data of the 
contained and upward going /i events, where #34 and 623 stand for, roughly speaking, the 
ratio of <-> v T versus <-> v s and the ratio of the contributions of sin 2 ( Am^^L/AE j 
versus sin 2 {Am^^L / AE) in the flavor oscillation probability, respectively, and 5\ is the 
only CP phase left in this scheme. Notice that the recent claim by the Superkamiokande 
group |27j that <-> u s is almost completely excluded is based on the hypothesis of two 
neutrino «-> v s oscillations with only one mass squared difference Am^ m , and their claim 
is completely consistent with the results in [Q, where the region 6*34 ~ ±7r/2 and #23 — 0, 
which would lead to pure <-> v s oscillations, is excluded. For generic mixing angles of #34 
and #23, the ansatz in |30| implies hybrid of Un «-> v T and z/„ <-> v s oscillations in general, and 



1 In the notation of [^J this parameter is given by c s = c\^c\^. We adopt different notation from 
[29] for the parametrization of the mixing angles, however, will use c s for |C4i| 2 + |t/s2| 2 - 
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one has to take into account the constraint of the CDHSW experiment [p8] also. Following 
30| , we take Am^ = Am1 SND = 0.3 eV 2 , which is consistent with the negative result of 



the disappearance experiment of CDHSW |28| for the entire region of the mixing angles 
obtained in |ID] . 

Combining the results of on solar neutrinos and [[JIJ on atmospheric neutrinos, we 
have evaluated the final ratio of the cosmic high-energy neutrino flux and the result is 
given in Figs. 4a and 4b, where the allowed region with the SMA (MSW), LMA (MSW) 
and VO solutions are shown separately. Almost the entire allowed region with the SMA 
(MSW) solution lies above the line F(u e ) = 1/3, and this is because the total normalization 
F{u e ) + F(is„) + F{y T ) = 1 - F(v s ) becomes less than 1 while F{u e ) = (1 - sin 2 26 l2 )F°{v e ) 
hardly changes due to smallness of = |0©|. On the other hand, for most of the region 
with LMA (MSW) and VO solutions, F{u e ) / F°{v e ) becomes smaller than 1 — F{y s ) and the 
region lies below the line F(u e ) = 1/3. This four neutrino scheme without BBN constraint 
gives us clearly a distinctive pattern for the final ratio of the cosmic high-energy neutrino 
flux, and observationally, if we have good precision then it may be even possible to distinguish 
the SMA (MSW) solution from the LMA (MSW) or VO solutions. The allowed regions in 
Figs. 4a and 4b are plotted for 5i = and 5± = tt/2. We observe that most of the allowed 
regions for 5i — and 5i = tt/2 overlap with each other and it implies that distinction 
between 5i = and 5i = tt/2 is difficult in this scheme of four neutrinos. 

We note in passing that the allowed region of the LOW solution to the solar netrino 
problem is contained in that of the LMA solution as far as sin 2 29 Q is concerned, and therefore 
the LOW solution gives us only a subset of the allowed region of the LMA solution in Figs. 
4a and 4b. 



III. NONSTANDARD RATIO OF THE HIGH-ENERGY NEUTRINO FLUX 

In section II, we have seen that the schemes of three neutrinos and of four neutrinos with 
BBN constraints give us the ratio F(i> e ): F{v^)\ F(u T ) ~ 1: 1: 1, irrespective of which solar 
solution is chosen. This is due to the fact that the intrinsic high-energy cosmic neutrino flux 
has the ratio F (u e ): F°(z/ M ): F°(u T ) = 1: 2: and the oscillation of atmospheric neutrinos 
is with maximal mixing while \U e ^\ 2 < 1. It has been pointed out || that the intrinsic 
flux of the cosmic high-energy neutrinos may not have the standard ratio F°{v e ): F°(^): 
F°(i/ T ) = 1: 2: 0, mainly because some of muons may lose their energy in a magnetic field. 
Here, we discuss in a model independent way the consequences of a generic scenario which 
is characterized by F (v e ): F°(z/ M ): F°(u T ) = A/3: 1 — A/3: 0, where A is a free parameter 
< A < 1 (the standard ratio is obtained for A = 1). 

We have plotted the allowed region in Figs. 5a and 5b for A = 1/2 and A = 0, respectively. 
We observe that the cosmic high-energy neutrino flux with the nonstandard ratio gives 
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relatively lower value of F(u e ). If relatively lower value of cosmic u e flux is observed in 
the future experiments, then it may imply that the possible oscillation scenario is the four 
neutrino scheme without BBN constraints and the intrinsic cosmic neutrino flux with the 
nonstandard ratio. Independent information on neutrino mixing parameters as well as on 
relevant astrophysical inputs may be needed here to arrive at a more definite conclusion. 

IV. RESULTS AND DISCUSSION 

In general, the final flux of high-energy cosmic neutrinos is expected to be almost equally 
distributed among the three (active) cosmic neutrino flavors because of vacuum flavor mix- 
ing/oscillations provided the astrophysical sources for these high-energy cosmic neutrinos 
are cosmologically distant, essentially irrespective of the neutrino flavors (three or four). 
Nevetheless, this may not be the case if the intrinsic high-energy cosmic neutrino flux ratios 
differ from the standard one, namely from F°(i/ e ): F ^): F°(u T ) = 1: 2: 0. In the exam- 
ples considered in this work for the nonstandard intrinsic high-energy cosmic neutrino flux, 
a relatively lower final flux for cosmic v e is obtained. The situation of nonstandard intrinsic 
cosmic neutrino flux may arise, for instance, if some of the muons lose their energy in the 
relatively intense magnetic field in the vicinity of the source. 

A simultaneous measurement of the three cosmic neutrino flavors may be useful to obtain 
information about a particular neutrino (mass and) mixing scheme depending on the relevant 
achievable resolutions for typical km 2 surface area neutrino telescopes. 

Irrespective of the numbers of neutrino flavors, in each of the neutrino (mass and) mixing 
scheme discussed in this work, the final flux of high-energy cosmic tau neutrinos is essen- 
tially comparable to that of non tau (active) neutrinos, even if it is intrinsically negligible. 
This may, at least in principle, be useful to constrain the relevant nonstandard particle 
physics/astro physics scenarios. 

In this work, we have considered the effects of vacuum neutrino flavor oscillations on 
high-energy cosmic neutrino flux in the context of three as well as four flavors. These 
oscillations result in an energy independent ratio, R a p = N a /Np (a ^ /3; a = (3 = u e , u T ) 
of the number of events detected for the neutrino flavors a and f3. It is so because the various 
neutrino flavor precession probabilities given in section II [see Eq. (5) and Eq. (17)] are 
essentially independent of neutrino energy. In the following paragraph, we briefly describe 
the prospects offered by the typical km 2 surface area under ice/water neutrino telescopes 
which are currently under construction/planning to possibly identify the cosmic neutrino 
flavor and hence to determine the ratio, R a p ||. 

We ignore the possible observational difference between cosmic neutrinos and anti neu- 
trinos for simplicity in the following discussion and assume that the flavor content in the 
cosmic neutrino flux is equally distributed because of vacuum flavor oscillations. Several 
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of the recent discussions suggest that the absorption of high-energy cosmic neutrino flux 
by earth is neutrino flavor dependent ||. The upward going electron and muon neutrino 
fluxes are significantly attenuated typically for Eq > 5 • 10 4 GeV, whereas the upward go- 
ing tau neutrinos with E > E may reach the detector with E < E because of the short 
life time of the associated tau lepton and may appear as a rather small pile up with fairly 
flat zenith angle dependence. For E > Eq, the upward going cosmic neutrino event rates 
range typically as: N u ~ (^(lO 1 ) whereas N Ut ~ (9(10°) in units of per year per steradian 
for typical km 2 surface area neutrino telescopes, if one uses the current upper high-energy 
cosmic neutrino flux limits ||. Let us note that for E > E , presently the high-enegy 
cosmic neutrino fluxes from AGNs dominate above the atmospheric neutrino background. 
For downward going high-enegy cosmic neutrino flux, the event rate ranges typically as: 
N Ve ~ O(10 1,5 ), N v ~ (9(10 2 ), whereas N Ut ~ (^(lO 1 ) in units of per year per steradian 
for the same high-energy cosmic neutrino fluxes. For E > Eq, the downward going cosmic 
tau neutrinos typically produce a two bang event topology such that the two bangs are con- 
nected by a /x-like track. The size of the second bang being on the average a factor of two 
larger than the first bang. The downward going electron neutrinos produce a single bang 
at these energies whereas the muon neutrinos typically produce a single shower alongwith a 
zipping /i-like track in km 2 surface area neutrino telescopes. Based on this rather distinct 
event topologies, cosmic neutrino flavor identification may be conceivable. The above order 
of magnitude estimates indicate that the typical km 2 surface area neutrino telescopes do 
offer some prospects for observations of high-energy cosmic neutrino flavor ratio, R a p, or at 
least may constrain it meaningfully. 
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Figures 



Fig.l The representation of the ratio of the final flux of (downward going) high-energy 
cosmic neutrinos on earth with a unit regular triangle. The x and y coordinates of the 
points are given by x = (2F T + F e ) / \/3, y = F e and the inside region of the triangle 
is given by < x < 2/ y/3, < y < 1. The point with an asterisk stands for the ratio 
without oscillations F(v e ): F{y^)\ F{y T ) = 1: 2: and is given by the coordinate 
(1/3^,1/3). 

Fig. 2 The ratio of the final flux of high-energy cosmic neutrinos in the far distance in the 
three neutrino scheme. The allowed region is the inside of each contour; (b) is an 
enlarged figure of (a). The allowed region lies near the mid point (l/\/3, 1/3). 

Fig. 3 The ratio of the final flux of high-energy cosmic neutrinos in the far distance in the 
four neutrino scheme with BBN constraints. The allowed region lies near the mid 
point (1/^3,1/3). 

Fig. 4 The ratio of the final flux of high-energy cosmic neutrinos in the far distance in the 
four neutrino scheme without BBN constraints; (b) is an enlarged figure of (a). The 
allowed region of the LOW solution is a subset of that of the LMA solution. 

Fig. 5 The ratio of the final flux of high-energy cosmic neutrinos in the far distance in a 
nonstandard scenario characterized by F°(u e ): F°(i/ (U ): F°(u T ) = A/3: 1 — A/3: 0. In 
Fig. 5(a), A = 1/2, the cases of N v = 3 with SMA and N v = 4 with BBN constraints 
have small region near the point F(u e ): F{y^: F(y T ) = 2: 5: 5, the region for the 
cases of N v = 3 with LMA and VO lie above this point, whereas most of the region for 
the case of N v = 4 without BBN constraints lie to the left of this point. As in Fig.4, 
the allowed region of the LOW solution is a subset of that of the LMA solution. In 
Fig. 5(b), A = 0, the cases of N v = 3 with SMA and N u = 4 with BBN constraints 
have small region near the point F(u e ): Fiy^): F(u T ) = 0: 1: 1, the region for the 
cases of N u = 3 with LMA and VO lie above this point, whereas most of the region 
for the case of N u = 4 without BBN constraints lie to the left of this point. 
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